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This  paper  addresses  the  dynamics  of  an  ultra-high  current  electron  ring  in  a  modified 
betatron  configuration.  Our  studies  include  analytical  and  numerical  results  for  both 
“cold”  and  “hot"  rings.  It  has  been  found  that  the  walls  surrounding  the  ring  and 
toroidal  effects  (hoop  stresses)  play  a  very  important  role  in  the  dynamics  of  the  ring, 
even  when  the  wall  conductivity  is  infinite.  For  finite  wall  conductivity,  the  diffusion  of 

(Continues) 


DO  .  JAN  n  1473  EOlTlON  OF  I  NOV  #S  If  OBSOLETE 

S/N  0  10  J-01«- «*01 


SECURITY  classification  OF  This  FAOE  (Ph«i  Dm ■  * MM) 


i 


Sltu«lT»  CLA$$t?lCATIOM  O*  Tmi«  HOC  f»—  0«1  <M 


I 

J^O.  AMTRACT  (Cmulmum0 

the  ring's  self  magnetic  field  profoundly  impacts  its  dynamics  and  in  general  the 
equilibrium  could  be  lost  if  means  are  not  provided  to  compensate  for  this  effect. 

In  addition,  it  has  been  found  that  the  Grad  B  drift  is  not  important,  except  when  the 
bounce  frequency  is  very  small. 

The  general  conclusion  of  these  studies  is  that  equilibrium  states  of  ultra-high  current 
rings  in  a  modified  betatron  configuration  exist  over  a  wide  range  of  parameters.  These 
states  are  accessible  and  within  the  realm  of  existing  pulsed  power  technology. 
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EQUILIBRIUM  OF  A  HIGH  CURRENT  ELECTRON  RING  IN  A 
MODIFIED  BETATRON  ACCELERATOR 

I.  Introduction 

Over  the  last  few  years,  several,  laboratories*”7  have  been  engaged  in  studies 
that  are  aimed  to  assess  the  feasibility  of  developing  ultra  high  current 
accelerators.  These  studies  are  mainly  motivated  by  the  potential  application  of 
high  energy,  high  current  electron  beams  in  nuclear  physics  research,  medical 
radiography  and  the  fusion  program. 

In  general,  the  effort  is  not  directed  toward  developing  novel  accelerating 
schemes  but  rather  in  modifying  or  improving  the  existing  accelerator  technology. 
Among  the  various  proposed  modifications,  the  addition  of  a  strong  toroidal  magnetic 
field  to  a  conventional  betatron®*^  has  attracted  considerable  attention^-7  and  the 
configuration  has  been  named  the  modified  betatron. 

In  this  paper  we  analyze  and  discuss  the  dynamics  of  a  high  current  electron 
ring  confined  in  a  modified  betatron  configuration.  When  the  intense  electron  ring 

Is  surrounded  by  a  finite  conductivity  wall,  its  dynamics  can  be  divided,  rather 

naturally,  into  three  distinct  phases:  The  pre-acceleration  phase  follows  the 
injection  and  trapping  of  the  electron  beam.  The  duration  of  this  phase  is  short, 
typically  a  few  psec  and  the  conducting  wall  surrounding  the  ring  can  be  treated  as 
perfect  conductor.  The  pre-acceleration  phase  is  followed  by  the  diffusion  phase, 
which  lasts  for  a  time  that  is  of  the  order  of  the  magnetic  field  diffusion  time. 
During  the  diffusion  phase  the  self  magnetic  field  of  the  beam  diffuses  out  of  the 
metal  torus.  As  a  result  the  electron  energy  is  reduced  and  the  hoop  stresses 

increase,  but  at  different  rates.  Therefore,  the  equilibrium  will  be  lost  if  means 

are  not  provided  to  balance  these  two  effects.  The  last  phase,  i.e. ,  the  main 
acceleration,  starts  after  the  self  magnetic  field  diffuses  out  of  the  torus  and  has 
a  duration  that  is  comparable  with  the  acceleration  time.  For  most  of  the  third 
Manuscript  submitted  August  17, 1982. 
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phase  Che  energy  of  Che  ring  has  increased  subscancially  and  Chus  Che  effecc  of  Che 
self  fields  is  appreciably  diminished. 

The  presenC  work  addresses  Che  dynamic?  of  high  currenc  ring  in  a  modified 

betatron  configuration.  It  includes  both  analyCical  and  computational  studies  for 

"cold”  as  well  as  "hot"  electron  rings.  The  main  conclusions  are:  First,  if  the 

energy  of  the  injected  electrons  is  not  exactly  equal  to  the  energy  corresponding  to 

the  equilibriun  orbit,  which  is  assumed  to  coincide  with  the  minor  axis  of  the 

torus,  the  center  of  the  orbit  is  displaced  from  the  center  of  the  minor  cross- 

section  of  the  torus.  The  displacement  is  proportional  to  the  energy  mismatch. 

This  imposes  very  stringent  constraints  on  the  injector.  Second,  during  the 

diffusion  of  the  self  magnetic  field  of  the  ring  out  of  the  chamber,  the  equilibrium 

can  be  lost  if  means  are  not  provided  to  balance  the  change  of  the  ring  equilibrium 

radius  that  are  due  to  the  reduction  of  the  relativistic  factor  y  and  the  increase 

o 

of  hoop  stresses.  Third,  the  modification  in  the  flux  rule  for  high  current  rings 

3 

can  be  ignored,  provided  that  v/yq  (v  *  Budker  parameter)  is  very  small.  Fourth, 
with  the  exception  of  the  expansion  of  the  minor  cross-section  of  the  ring,  finite 
emittance  does  not  have  any  other  noticable  effect  on  the  equilibrium,  and  fifth, 
considerably  higher  axial  energy  spread  can  be  tolerated  in  high  than  in  low  current 
rings. 

The  general  and  most  important  conclusion  of  these  studies  is  that  equilibrium 
states  of  high  current  rings  in  a  modified  betatron  configuration  exist  over  a  wide 
range  of  parameters.  These  equilibria  are  realistic  and  accessible  with  existing 
technology. 
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II.  Orbic  Stability  of  a  Cold  Ring  Without  Surrounding  Walls  and  Toroidal 


Corrections 

The  stability  of  a  single  particle  orbit,  when  thermal  effects,  wall  effects 
and  toroidal  corrections  are  neglected,  has  been  considered  previously.^  Using  the 
coordinate  system  shown  in  Fig.  1  and  assuming  that  the  external  field  components 
vary  as 


B  (r ,t)  -  B  ( t ) 
z  oz 

{l  -  n  (r  -  ro)/rQ}, 

(la) 

Br(r,t)  -  - 

Boz(t)  °*/ro» 

(lb) 

E0  (r-c)  -  -  TZ  'l 

rt 

3B 

r'dr'  (r\  t) , 

(lc) 

o 


and 


Be(r)  "  Boe  ^  "  (r  “  W’ 


(Id) 


it  has  been  determined  that  for  Bfl  »  B  and  v/y  <  <  1  the  orbits  are  stable, 

“  z  o 

provided  that  the  following  condition  is  satisfied 

"s  i  "W28OI>2-  <2> 

In  the  above  equations  B  (r,t)  is  the  axial  and  B  (r,t)  is  the  radial  component  of 
the  betatron  field,  B0(r)  is  the  toroidal  magnetic  field,  E0(r,t)  is  the  induced 


electric  field,  a  Is  the  external  field  index,  rQ,zo  are  the  coordinates  of  the 

center  of  the  electron  ring  and  n  is  the  self-field  index.  The  self-field  index  is 

s 

defined  as 

",  •  “t,2'‘2Vo*2>  '  (c/ao,tb>2-  <3) 

2  2 

where  {*  4ire  nQ/m}  is  the  beam  plasma  frequency  squared,  yq  is  the  relativistic 
factor,  Q  {*eBoz/mc}  is  the  cyclotron  frequency  corresponding  to  the  axial 
component  of  the  betatron  field,  r.  is  the  minor  radius  of  the  beam  and  v  is 

D 

Budker's  parameter,  i.e.,  the  product  of  the  number  of  electrons  per  unit  length 
times  the  electron  classical  radius. 

The  electron  beam  current  that  can  be  stably  confined  in  a  modified  betatron 
can  be  obtained  by  substituting  Eq.  (3)  into  Eq.  (2)  an \  is 

‘.b  S  2-‘  <rb/ro>2  *.V  <Boi/V2  <“>•  (4) 

where  0  ■  v  /c.  When  y  »  1,  then  r  0  /y  ■  c  and  Eq.  (4)  becomes 

o  o  o  o  oz  o 

_7  2  2 

1  .  <  7.22  x  10  r.  y  B  Q  {ka},  with  r,  in  cm  and  B  Q  in  Gauss, 
mo  —  oooo  0  oo 

In  addition,  it  has  been  shown  under  the  same  conditions  but  with  »  0 
(conventional  betatron) ,  the  beam  is  stable  provided 

n  <  1/2,  conventional  betatron.  (5) 

3  <b8  -  °> 
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The  electron  beam  current  that  can  be  stably  confined  in  a  conventional 
betatron  can  be  obtained  by  substituting  Eq.  (3)  into  Eq.  (3)  and  is 

Icb  <  (wW  w 

or 


I  .  <  14.44  x  10"7  r  2y  B  2  IkaI, 

Cb  -  K  n  r\T  i  * 


b  o  oz 


when  y  »  1.  The  ratio  of  I  .  to  1  .  obtained  from  Eqs.  (4)  and  (6)  is 
o  no  cd 


rr- 1/2  (B0e/Boz)2* 

cb 


The  above  relation  indicates  that  for  B  .  »  B  ,  the  electron  beam  current  that  can 

o  8  oz 

be  stably  confined  in  a  modified  betatron  substantially  exceeds  the  current  that  can 
be  confined  in  a  conventional  betatron. 

The  modified  betatron  stability  condition  given  in  Eq.  (2)  can  be  easily 
obtained  from  the  well  known  confined  equilibrium  condition^® 

2(  tu^2/^2)  <  1,  (non-relativistic)  (7) 


where  <0^  is  the  beam  plasma  frequency  and  is  the  cyclotron  frequency.  For 
relativistic  energies  and  taking  into  account  the  self  magnetic  field  of  the  beam 
becomes 


2 

% 


_1_ 

,  2’ 

% 

self  field  correction 


(8) 


6 


and 


c 


(9) 


Substituting  Eqs.  (8)  and  (9)  into  Eq.  (7)  we  obtain  the  stability  condition  of  Eq. 

(2).  The  equality  sign  in  Eq.  (2)  gives  the  maximum  electron  density  that  ca~  be 

supported  at  a  specific  value  of  B  ./ B  and  the  corresponding  equilibria  i  nown 

o  o  oz 

as  Brillouin  flow.*® 


III.  Mall  Effects  on  the  Macroscopic  Motion  of  a  Cold  Beam 

In  this  section  we  analyze  the  effects  of  surrounding  walls  on  the  motion  of 
the  center  of  the  beam.  In  the  first  sub-section,  it  is  assumed  that  the  perfect 
conductor  that  surrounds  the  beam  is  a  straight  cylindrical  pipe  of  circular  cross- 
section  and  thus  toroidal  effects  (hoop  stresses)  are  omitted.  These  effects  are 
included  in  the  second  sub-section. 


a.  Without  toroidal  corrections 

Consider  a  pencil-like  electron  beam  inside  a  straight,  perfectly  conducting 
cylindrical  pipe  of  circular  cross-section  as  shown  in  Fig.  2.  The  center  of  the 
beam  is  located  at  a  distance  A r,  tz  from  the  center  of  the  minor  cross-section  of 
the  pipe.  As  a  result  of  the  induced  charges  on  the  wall,  the  center  of  the  beam 
will  experience  a  radial,  outward  directed  force,  which  for  small  displacements, 
i.e.,Ar,  Az  «  a  is  given  by 


F£  -  2ire2nQ(rb/a)2  {Arer  +  Azez}, 


(10) 
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Figure  2. 


Wall  (images)  forces  acting  on  a  pencil-like  electron  beam 
situated  inside  a  perfectly  conducting  cylindrical  pipe. 


where  a  is  the  cylinder  radius  and  nQ  the  uniform  beam  density. 

Similarly,  as  a  result  of  the  induced  current  on  the  wall,  the  center  of  the 

beam  will  experience  a  radial  force  that  is  directed  toward  the  opposite  direction 

than  Fe  and  is  given  by 
E 


f»  ■  '  8o2  FV 


Using  the  external  fields  of  Eq.  (1)  and  the  induced  fields  of  Eqs.  (10)  and  (11), 
the  equations  describing  the  temporal  linear  evolution  of  the  beam's  center,  for 
time  independent  applied  fields,  are: 


•*  ~  2 
Ar  +  u  Ar 
r 


a  a  0  <5P0> 

-21  Az*  +  r-2£)  - 2- 

Y  vy  '  y  mr 

o  o  o  o 


-  2 

Az  +  a:  Az 

z 


2  f2 

where  m.2  -  '  n  *  ns  ~r}  *  V  "  (noz/Yo)2(n  “  ns  ”T)  * 

a  a 

«Y  0  <  5P  g> 

Q  „  •  eB  Jmc,  a  -  eB  /me,  -  *  -  and  6Pa  is  the  difference  between  the 

0  0  off  ’  oz  oz  ’  Y  Y  mr  c  9 

o  oo 

canonical  angular  momentum  of  an  electron  at  (r,z)  and  its  corresponding  value  at 
the  equilibrium  orbit  (rQ,o).  The  average  is  over  initial  coordinates  and 
velocities.  Equations  (12)  and  (13)  do  not  include  the  self  electric  and  self 
magnetic  fields,  because  both  these  fields  are  zero  at  the  center  of  a  straight 
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beam.  In  addition,  the  non-linear  terms - Am  and - —  have  been  omitted 

r0  Yo  r0  Yo 

from  Eqs.  (12)  and  (13).  These  two  terms  have  their  origin  in  the  gradient  of  the 
toroidal  magnetic  field  and  are  considered  in  Section  V.  In  general  these  terms  are 
not  significant  except  in  the  limit  ♦  o. 

In  Eq.  (12),  ■  80<5P  g>/mrQc  indicates  the  energy  mismatch,  i.e.  ,  the 

difference  between  the  energy  of  the  reference  electron  (moving  along  the  axis  of 
the  beam)  and  the  energy  required  for  the  same  electron  to  move  on  the  equilibrium 
orbit  (rQ,o).  The  solution  of  Eqs.  (12)  and  (13),  for  time  independent  fields,  is 


Ar  ■ 


noz<6pe> 
2-  2 

y  a)  r  m 
o  r  o 


,4  ,-  2  2.1/2  .  „ 

+  j*i  cj  (“2  '  “j  >  .‘y* 


(14) 


and 


Az 


A 


;i  (uj 


'  2.1/2 

-  »r  ) 


(15) 


where  c 


j 


is  a  constant  and 


-  1/2  {w.2  +  a,2  +  (no0/Yo)  +  [(wrZ  +  + 


-2.-2 


(16) 


5«2il/2^ 


The  first  term  on  the  RHS  of  Eq.  (14)  gives  the  displacement  of  the  center  of 
the  orbit  from  the  center  of  the  surrounding  cylindrical  pipe  and  can  be  written  as 


Ar 

o 


<*V/Y0 _ 

(VV  (l  “ n  •  Vb2/a2)mro 


(1?) 
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Equation  (16)  has  four  roots.  TVo  of  them  are  fast  (plus  sign)  and  two  are  slow 


(minus  sign). 


When  BQ 
^  2 


0)  u> 

r  z 

^o6/ro 


»  B  ,  the  slow  mode  becomes 
z 


(18) 


and  the  fast  mode  is  *  (Q  Jy  ).  Equation  (18)  is  plotted  in  Fig.  3  for  two  values 

O  0  o 

2 


n  rw 

s  b 


of  external  field  index.  For  n  >  1/2,  is  negative  when  1  -  n  <  — - —  <  n  and 

a 

2 


n  r. 
s  b 


for  n  <  1/2,  Wg  is  negative  when  n  <  — ^ —  <  1  -  n.  Negative  values 


of  <d  indicate  that  the  beam  motion  is  unstable  and  the  orbit  in  the  r,z  plane  is 

D 


2  2  -3 

open.  Since  the  parameter  ngrb  / a  scales  as  yq  ,  during 


2,2 


acceleration  ogr^  /a  decreases  rapidly.  Therefore,  in  order  to  avoid  the 


instability,  it  is  necessary  that  before  the  commencement  of  the  acceleration  the 


2/2 


parameter  ngr^  /a  is  less  than  1  -  n  when  n  >  1/2  and  less  than  n  when  n  <  1/2. 


This  implies  that  the  injected  beam  current  should  be  limited  to 


I  <  8.5  (1  -  n)  8  y  3a2/r  2  {ka}  for  n  >  1/2 
—  0  0  0 


and 


I  <  8.5  nB  y  3a2/r  2  {ka}  for  n  <  1/2. 
—  o  o  o 


It  is  apparent  that  these  two  constraints  are  rather  lenient. 


The  orbit  of  the  ring's  center  is  described  by 


Ar(t)  ■  [Ar(o)  -  Ar  ]  cos  ui_t  +  sin  u  t  +  Ar  , 

O  o  (4^  B  O 


Az(  t)  -  Ar(o)  cos  (<ut)  +  Az--°--  sin  (u>  t), 

B  “b  B 


where  the  initial  velocities  and  displacements  are  related  by 


Ar(o)  -  -  [ujz2/CnoQ/ yq)  ]  Az(o), 


,  n  <5Pq> 

Az(o)  -  [Sr2/(no9/YQ)]  Ar(o)  + 

00  o  o 


and 


<5Pfl>  a  rnarJva~  rn 

t)  O  O  8  0  oz  o 


The  predictions  of  Eq.  (18)  are  in  good  agreement  with  the  results  of  computer 
simulation  shown  in  Fig.  A.  This  figure  shows  four  snap-shots  of  the  beam  in  the 
r,z  plane.  At  t  *  0  the  circular  electron  beam  is  injected  near  the  center  of  the 
pipe.  The  values  of  the  various  parameters  are  listed  in  Table  I.  The  electron 
beam  current  is  kept  low  (1  KA)  in  order  to  minimize  the  toroidal  effects.  As  can 
be  seen  from  Fig.  A,  the  center  of  the  electron  beam  describes  a  circle  in  the  r,z 
plane  with  a  period  of  about  188  nsec.  For  the  same  parameters  Eq.  (18)  predicts  a 
period  ig  ■  2ir/aj ^  -  180  nsec.  These  numerical  results  uce  discussed  further  in  the 
next  sub-section. 

The  displacement  of  the  orbit's  center  because  of  the  energy  mismatch  [Eq. 

(17)],  imposes  very  stringent  constraints  on  the  injector.  This  becomes  apparent 

2 

when  we  consider  some  limiting  cases.  For  example,  when  n  =  4  and  <<  1,  Eq. 

(17)  is  reduced  to 


13 


(cm) 


43.6 


R  (cm) 


Figure  4. 


Snap-shots  of  the  electron  ring  minor  cross-section  at  foul 
different  times.  The  values  of  the  various  parameters  are 
in  Table  I.  The  center  of  the  ring's  minor  cross-section 
describes  a  circle. 
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Table  I. 


Parameters _ 

Beam  Energy  (MeV)  3 

Beam  Current  (kA)  1 

Beam  Minor  Radius  (cm)  1 

Beam  Major  Radius  (cm)  50 

Torus  minor  radius  (cm)  6.4 

Vertical  magnetic  field  (G)  240 

Toroidal  magnetic  field  (kG)  2 

External  field  index  0.47 

Self  field  index  0.88 


15 


(19) 


2(.Sv/v). 

0  o 


Equation  (19)  predicts  that  for  a  major  radius  vq  *  100  cm,  the  ratio  6  yq  /  YQ  should 
be  less  than  1%  in  order  the  displacement  of  the  orbit  to  be  less  than  2  cm.  The 
condition  6yq/ yq  <  1%  requires  that  the  uncertainty  in  energy  should  be  less  than  35 
KeV,  when  the  energy  of  the  injected  beam  is  3  MeV.  Although  such  a  small  energy 
uncertainty  can  be  attained  with  sophisticated  injectors,  there  are  other  factors, 
such  as  space  charge  and  inductive  effects,  which  contribute  substantially  to  the 
uncertainty  of  energy.  Although  the  displacement  cannot  be  eliminated,  its  negative 
impact  can  be  alleviated  by  elongating  the  minor  cross-section  of  the  torus  along 
the  r-direction. 


b.  Toroidal  Corrections^’11 


The  cause  of  these  effects  is  the  finite  curvature  of  the  electron  beam 
orbit.  For  relative  large  aspect  ratio  r^r^  »  1  beams,  the  toroidal  effects 
become  important  when  v/ yq  exceeds  a  few  percent.  Previous  work1^  on  toroidal 
effects  was  limited  to  "bare"  rings,  i.e.,  without  surrounding  conducting  walls.  In 
this  sub-section  we  briefly  analyze  the  toroidal  corrections  for  the  more  realistic 
geometry  shewn  in  Fig.  5,  i.e.,  including  the  effect  of  conducting  wall  around  the 
beam.  Since  different  physics  issues  are  involved  in  the  case  of  a  "bare"  and 
"shielded"  rings,  it  is  appropriate  to  start  our  discussion  with  a  bare  ring. 

Consider  on  intense  electron  ring  as  shown  in  Fig.  5,  but  without  the  toroidal 
chamber.  At  the  inner  edge  of  the  RHS  cross-section  of  the  ring  the  self-magnetic 
field  is  greater  than  that  of  a  straight  beam  with  the  same  parameters,  because  of 
the  contribution  from  the  LHS  of  the  ring.  At  the  outer  edge  of  the  cross-section, 
the  self  magnetic  field  is  reduced  because  the  contribution  from  the  LHS  of  the  ring 
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Figure  S.  An  electron  ring  inside  a  conducting  toroidal  chamber.  Because  of 
the  toroidal  effects,  neither  the  self-electric  not  the  self 
magnetic  field  is  zero  at  the  geometric  center  of  the  ring. 


has  opposite  polarity  than  the  local  field.  Thus,  the  total  self  magnetic  field  is 
shifted  upwards,  as  shown  schematically  in  Fig.  5.  In  contrast,  the  self  electric 
field  decreases  at  the  inner  radius  and  its  magnitude  increases  at  the  outer  radius 
of  the  ring.  Thus  the  self-electric  field  shifts  downwards.  As  a  result  of  these 
shifts  neither  the  magnetic  nor  the  electric  field  are  zero  at  the  geometric  center 
of  the  ring. 

Uhen  the  electron  ring  is  surrounded  by  a  perfect  conductor,  the  shift  in  the 
fields  discussed  above  is  reduced  but  additional  field  components  appear  as  a  result 
of  induced  charge  and  current  on  the  conducting  wall.  For  small  ring  displacements 
from  the  center  of  the  minor  cross  section  of  the  torus,  the  induced  fields  vary 
linearly  with  the  displacement  and  are  identical  to  those  given  in  Eqs.  (10)  and 
(11),  for  a  cylindrical  pipe.  Quantitative  extension  of  the  previous  qualitative 
considerations  has  shown  that  the  fields  at  the  center  of  an  uniform  density 
electron  ring  inside  a  perfectly  conducting  toroidal  chamber  of  circular  cross- 
section  are^’^ 


~  ind 


•2ir|e|n  r 
o  o 


r,rb  Ar  1  rb  a  .  ~  rb  Az 

l(T  r  +  IT  talT)  er  +~- 

a  o  r  b  a  o 


and 


-2x|  e|n  8  r  [4-  ±  e  - 

o  o  o  L  2  r  r 


«*♦!*"  f-l>  «.]• 

a  o  r  b 


where  n  is  the  ambient  density,  8  ■  v  /c,  v  is  the  azimuthal  velocity  defined  by 

O  OO  O  ' 

r  Jl  fy 

_ o  oz  o _  f?7 ■> 

Vo  *  1  +  2(v/y  )(1  +  in  a/rw)  ’  1 

0  D 
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and  Che  displacement  Ar ,  ts.  of  Che  ring  from  Che  center  of  the  torus  has  been 
assumed  to  be  much  less  than  a.  In  addition,  the  fields  given  by  Eqs.  (20)  and  (21) 
have  been  derived  under  the  assumption  that  the  angular  frequency  of  the  electrons 
is  constant  and  therefore  the  electron  current  density  varies  proportionally  to  r. 
Using  Eqs.  (20)-(22),  it  can  be  shown  that  the  center  of  the  beam  is  described 

^  2  2 

by  Eqs*  (12)  and  (13)  with  and  replaced  by 


2/  2, 

-  ngrb  /a  ), 


2/  2. 
Vb  /a  >• 


(23) 


(24) 


and 


«p9  -  {*, 


where 


a  -  (1 
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€ 


[i  +^a  +  «n|-)] 
7o  b 


-1 


(25) 


and 


* 

n  -  n£. 


The  bounce  frequency  can  be  found  by  substituting  Eqs.  (23)  and  (24)  into  Eq. 


(18)  and  is 


(26) 


When  n*  ■  a/2,  for  the  reasons  stated  In  the  paragraphs  following  Eq.  (18),  the 
orbits  are  closed  (stable)  as  yq  increases,  provided 

n  r.2/a2  <  a/2 

S  D 

or,  for  y  »  1 
o 

I  <  4-25  (f-)2  (1  -  ~  In  f-)  Y03  [ka].  (27) 

b  o  b 

The  limiting  current  given  by  Eq.  (27)  is  considerably  lower  than  that  of  Eq. 

(4)  when  B  n  »  B  .  In  addition,  in  contrast  to  the  current  given  by  Eq.  (4),  the 
o  o  oz 

current  of  Eq.  (27)  is  independent  of  the  toroidal  magnetic  field  Bq0. 

At  this  point,  it  is  appropriate  to  return  and  discuss  further  the  computer 

simulation  results  of  Fig.  4.  For  the  l  KA,  3  MeV  beam  the  v/ yq  -  0.0084.  Even  at 

this  small  value  of  v/y  toroidal  effects  are  noticeable.  The  value  of  £  computed 

o 

from  Eq.  (25)  using  the  above  value  of  v/ yq  and  6.4  for  the  ratio  a/r^  is  0.95. 

According  to  Eqs.  (23)  and  (24)  the  center  of  the  beam  will  describe  a  circle 

when  a  -  n  -  n  or  for  n  *  :kl  -  v/ Y  to  a/r,  )£  ■  0.467,  which  is  in  excellent 

l  ob 

agreement  with  the  simulation.  In  addition  Eq.  (26)  predicts  a  period  of  188  nsec, 
which  also  is  in  excellent  agreement  with  the  simulation. 

The  most  striking  manifestation  of  toroidal  effects  is  in  the  value  of  betatron 
magnetic  field  required  to  confine  the  rotating  beam  at  a  specific  radius.  When  the 
axis  of  the  beam  lies  along  the  axis  of  the  torus,  i.e.,  when  dr  ■  bz  ■  0,  it  can  be 
shown  from  Eqs.  (20)  and  (21)  that  the  external  magnetic  field  required  for  the  beam 
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Co  rocace  with  a  radius  r  is 

o 

B  "  I1  +  2u/Yrt  (1  +  In  a/r.)}, 
oz  o  o  b  ‘ 

where  Bq  is  the  magnetic  field  necessary  for  a  single  particle  of  the  same  energy  to 

rotate  with  a  radius  r  . 

o 

The  above  expression  for  the  magnetic  field  is  based  on  the  assumption  that  all 
the  electrons  rotate  with  a  constant  angular  frequency,  i.e.,  the  current  density 
increases  linearly  across  the  beam.  If  the  current  density  is  constant  across  the 
beam,  the  above  expression  is  slightly  modified  and  becomes 


B  -  Brt  {l  +  2v/y  (0.5  +  In  a/r.)}. 
oz  0  1  o  b  ' 


For  a  10  KA,  2  MeV  uniform  current  density  beam  with  a  ratio  a/r^  “  6.4,  the 
correction  is  55Z,  i.e.,  B  /B  *  1.55.  This  effect  is  demonstrated  clearly  by  the 
results  of  computer  simulation  shown  in  Fig.  6.  The  three  snap-shots  of  the 
electron  ring  minor  cross  section  in  a  modified  betatron  field  correspond 
to  t  *  0,  20  and  40  nsec.  For  all  practical  purposes,  the  minor  radius  of  the  beam 
remains  constant.  The  external  betatron  magnetic  field  is  127  G,  i.e.,  ~  50%  higher 
than  the  single  particle  field.  The  remaining  parameters  are  summarized  in  Table 
II. 

Another  consequence  of  the  toroidal  effects  is  the  increase  in  the  radial 

displacement  of  the  orbit  for  fixed  energy  mismatch  and  major  radius.  For 

n  *  a/2  and  n  r.^/a^  «  1,  Eq.  (19)  becomes  dr  /r  *  2(6y  / Y  )  (1  +  —  In  — ) .  In 
sb  ooooYtt 

2  2  2  o  b 

addition,  the  w  vs.  nr  / a  curves  of  Fig.  3  are  shifted  to  the  left  and  thus  the 
B  s  b 
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Figure  6.  Three  snap-shots  of  the  electron  ring  minor  cross  section  in  a 

modified  betatron  field.  The  external  magnetic  field  required  to 
confine  the  ring  is  502  higher  than  the  single  particle  field. 

The  values  of  the  various  parameters  in  this  run  are  listed  in 
Table  II. 


Table  II. 

Parameters  of  the  uniform  current  density  ring 


Beam  Energy  (MeV) 

Beam  Current  (kA) 

Beam  minor  radius  (cm) 

Beam  major  radius  (cm) 

Torus  minor  radius  (cm) 
Vertical  magnetic  field  (G) 
Toroidal  magnetic  field  (kG) 
External  field  index 


2 

10 

1 

100 

6.4 


127 


1.7 

0.25 
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maximum  permissible  current  that  can  be  accelerated  is  reduced,  as  may  be  seen  from 
Eq.  (27). 

IV.  Acceleration 

After  injection  and  trapping  the  beam  is  accelerated  to  high  energy  by  the 
inductive  electric  field  generated  by  the  time  varying  betatron  magnetic  field.  As 
shown  in  Fig.  7,  the  acceleration  may  be  divided  into  three  phases 

1.  Pre-acceleration 

2.  Diffusion 

3.  Main  acceleration 

The  pre-acceleration  phase  occurs  for  times  much  shorter  than  the  field 

4ir  r 

diffusion  time,  i.e.,  t  «  t  -  •  a  Sain  (— )  ,  where  a  is  the  conductivity,  6  is 

U  L  3 

c 

the  thickness,  a  is  the  minor  radius  and  rQ  the  major  radius  of  the  torus.  During 
the  pre-acceleration  phase  the  self  magnetic  field  of  the  beam  does  not  have  time  to 
diffuse  out  of  the  conducting  torus.  In  the  example  shown  in  Fig.  7  the  total 
acceleration  time  was  chosen  I  msec  and  the  diffusion  time  10  ysec.  The  ratio  of 
the  temporal  extent  of  each  phase  ^o  the  total  acceleration  time  t/t  is  also  given 

dL 

in  the  figure.  It  is  apparent  that  the  beam  spends  most  of  its  time  in  the  main 
acceleration  phase. 

During  the  pre-acceleration  phase  the  metal  wall  surrounding  the  beam  can  be 
treated  as  a  perfect  conductor.  As  a  result,  for  small  displacement  of  the  beam, 
the  self  flux  linking  the  axis  of  the  electron  ring  remain  constant  to  lowest 
order.  This  point  is  discussed  further  later. 

The  condition  for  the  major  radius  of  the  electron  ring  to  remain  constant  can 
be  derived  from  the  conservation  of  canonical  angular  momentum 
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6Xt 

where  v  is  the  toroidal  velocity,  is  the  vector  potential  that  describes  the 

o  o  o 

external  magnetic  (betatron)  field,  Ag  is  the  vector  pote  tial  that  describes  the 
self  magnetic  field  of  the  beam  and  r  is  the  instantaneous  radius  of  the  beam. 

The  self  flux  linking  the  ax.s  of  the  ring  $s  is 


4S  -  2xrAgS , 


(29) 


and  to  lowest  order  is  independent  of  the  beam  displacement  and  thus  of  time  See 
Eq.  (35)  .  Substituting  Eq.  (29)  into  Eq .  (28)  and  using  the  expression  for 


rfl  ly 
oz  o 

'o9  ”  1  +  2  T v/ y  )(l  +  In  a7r.  )  * 

o  b 


Eq.  (28)  predicts  that  r  remains  constant,  provided 
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(50) 


In  equation  (30),  BQz(t)  is  local  magnetic  field,  v  is  the  Budker' s  parameter, 

gXt 

a  is  the  minor  radius  of  the  torus,  is  the  beam  radius  and  <B  >  is  the  average 


applied  magnetic  field. 
Using  the  equation 


together  with 


I_  i_  <Rext> 
2c  at  <B  >’ 


(31b) 


equation  (30)  becomes 


3B  (t) 
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at 


1  U 

2  ]1 
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(1  +  in 
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r,  ; '  at 
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<BeXt>. 


(32) 


1  3 

A  similar  condition  has  been  derived  previously  . 

Equation  (32)  is  the  condition  that  must  be  satisfied  in  order  for  the  radius 
of  the  accelerated  electron  ring  to  remain  constant.  For  low  v/yq  beams,  Eq.  (32) 
is  reduced  to  the  well  known  flux  rule,  i.e., 


3B 


oz 


(t) 


at 


1  JL  cnext 

2  at  <B 


>. 
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The  correction  term  — j  (1  +  in  — )  in  Eq.  (32)  is  very  sensitive  to  the  beam 

Y  J  rb 

energy.  For  a  10  kA,  3  Mev  beam  injected  into  a  10  cm  minor  radius  torus  with 

r,  ■  1  cm,  the  correction  is  onlv  1%  and  therefore  can  be  neglected.  However,  when 
b 

the  energy  of  the  same  beam  is  reduced  to  0.5  MeV  the  correction  is  48%,  i.e.,  very 

substantial.  The  design  of  the  accelerator  is  simplified  considerably  by  choosing 

the  beam  parameters  such  that  the  correction  terra  is  negligible. 

The  instantaneous  value  of  y  (t)  may  be  determined  from  Eqs.  (30)  and  (31)  and 

o 

for  y  (t)  >>  1  is 
o 

Y  (t)  [l  +  2 v/ y  (0) ( 1  +  in  a/r,)]  B  (t) 
o  o _ b  oz _ 

(0)  [1  +  2 v/ YQ( t)  (1  +  in  a/rfe)  I  g 

0  oz 

In  contrast  to  the  lew  current  beams,  in  which  y  (t)/B  (t)  remains  approximated 

o  oz 
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constant  during  acceleration,  for  high  current  beams  with  toroidal  corrections  the 
quantity  that  remains  constant  during  acceleration  is 


~ — 7~y  (l  +  2v/yq(c)  (1  +  in  a/rfe)}  *  constant.  (33) 

oz '  0 


As  YQ(t)  increases  with  time,  the  correction  term  decreases  and  Eq.  (33)  is 

reduced  y  (t)/B  (t)  -  const, 

o  oz 

Now,  we  return  to  discuss  the  self  flux  linking  the  axis  of  the  ring.  In  the 
notation  of  the  toroidal  geometry  shown  in  Fig.  8a,  the  magnetic  stream  function 

S  11 

*  rA.  for  p  <  r,  is  given,  to  lowest  order,  by 

0  D 


1>/Q  *  2  in  a/rfe  +  (1  -  p2/rfe2)  -  2  p/a  j  cos  (i|>  -  <k)  +  (  in  ^ - + 

3  Kb  rb 


(34) 


2  2 

where  Q  ■  -cec  n  xr.  ftR.  /c, 

O  D  D 


vQ/r,A 


Ar2  + 


Az 


is  the  uniform  density  and 


it  has  been  assumed  that  a  «  r  ,  r.  «  a  and  Ar,  Az  <<  a. 

o  b  ’ 

The  self-flux  through  the  axis  of  ring  is 


$S  *  2x«K  p  ■  0)  . 


(35) 


Substituting  Eq.  (34)  into  Eq.  (35),  we  get 


<fs  -  2*q  (1  +  2  in  a/rb),  (36) 


which  does  not  depend  on  the  displacement  Ar,  Az  of  the  beam.  This  is  a  rather 
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unusal  result  and  deserves  further  discussion. 

To  gain  some  insight  into  the  problem,  we  have  computed  the  flux  linking  a 
horizontal  surface  s  extended  from  the  axis  of  a  straight  beam  to  the  inner  wall  of 
a  perfectly  conducting  cylinder  of  circular  cross-section  as  shown  in  Fig.  8b.  The 
beam  displacement  is  arbitrary  but  the  ratio  r^/a  «  1.  It  is  straight  forward  to 
show  that  the  flux  through  such  a  surface  of  length  l  is 

2  2 

*S  -  (2 11/ c)  {to  (a/r  )  -  ln(l  -  *  *  *  )}  (37) 

a 

where  I  is  the  beam  current. 

For  (Ar  +  Az^;  /a  «  1,  Eq .  (37)  indicates  that  the  flux  <fs  has  a  quadratic 
dependence  on  the  beam  displacement.  Therefore,  since  only  the  linear  terms  on  beam 
displacement  were  kept  in  the  derivation  of  Eq.  (34),  it  is  not  surprising  that  the 
flux  given  by  Eq.  (36)  is  independent  of  the  displacement  of  the  beam. 

When  the  time  approaches  the  magnetic  field  penetration  time  t^,  the  self 
magnetic  field  of  the  electron  beam  starts  to  diffuse  out  of  the  finite  conductivity 
metal  torus.  Using  the  geometry  shown  in  Fig.  9,  it  is  shown  in  Appendix  A,  that 
for  a  very  thin  conductor,  i.e. ,  a  *  b,  to  lowest  order,  the  self  magnetic  field  of 
the  beam  at  time  t  is 


-t/tT 


B<t  (r,t)  “cr^1-6  °  3 »  r  2.  b 


B*  (r*C)  "  «  *  rb  <  r  <  a> 


(38) 


where 
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Figure  3b.  A  displaced  straight  electron  beam  inside  a  perfectly  conducting 

pipe.  The  image  current  is  located  at  a  distance 
2  2  2 

a  /(Ar  +  Az  )  from  the  center  of  the  pipe. 
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These  equations  are  based  on  the  assumption 

that  6/a  «  1,  ((6/b)in  — «  1  and  the  inductive  electric  field  varies 

b  3Vb»c)  r 

as  Ez(r,t)  -  —  — -  in  (“) .  i-e.,  goes  to  zero  at  r  *  rQ  (pseudo-toroidal 


geometry) . 

The  self  magnetic  flux  in  the  three  regions  is  given  by 


4xr 

- -  Iin(a/r.) 

C  0 


r  <  a 


a  <  r  <  b 


4irr  -t/t_  r 

— 2l[l-e 


r  >  b 


and  it  diffuses  out  at  the  rate 


4irr  I 
o 

S' 


in  (— ) . 


(39) 


The  inductive  electric  field  generated  by  the  changing  flux  given  in  Eq.  (39) 

acts  to  slow  down  the  beam.  In  addition,  for  a  constant  current  density  ring,  the 

8r 

hoop  stresses  increase  by  the  ratio  (2  +  in  - )/(l  +  in  — ),  and  the  induced 

erb  rb 

magnetic  field  components  go  to  zero  at  the  end  of  the  diffusion.  However,  the 
electric  field  components  remain  the  same.  As  a  result,  for  a  electron  ring  that  is 
situated  along  the  minor  axis  of  the  torus,  the  decrease  in  its  equilibrium  radius 
associated  with  the  reduction  of  is  greater  than  the  corresponding  increase  of 
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Che  equilibrium  radius  associated  with  the  enhanced  hoop  stresses  and  thus 
equilibrium  can  be  lost.  This  difficulty  can  be  avoided  by  placing  a  set  of 
external  conductors  along  the  minor  cross-section  of  the  torus  having  a  poloidal 
distribution  that  closely  resembles  the  distribution  of  wall  currents  in  a  perfect 
conductor,  i.e., 


I 

w 


*2 w  "  zfr 


where  ^  is  the  wall  current  per  unit  length  and  1^  the  ring  current. 

This  compensation  is  satisfactory  even  when  the  beam  is  displaced  off  center 
and  rotates  around  the  equilibrium  position  with  provided  that  »  1*  The 

reason  is  that  the  correction  term  in  the  fields  at  the  center  of  the  ring,14  when 
the  skin  depth  is  much  greated  than  the  thickness  of  the  conductor,  is 


—  t  /  T 


-  2  e 


UL  T 
B  D 


sin  u)  t, 


and  therefore,  it  can  be  neglected.  Similarly,  when  the  skin  depth  is  much  smaller 
than  the  thickness,  the  correction  term  is  also  small  and  is  given  by 


-t/ T, 


2  e 


T 

B  D 


s  inui^t. 


In  addition,  we  have  shown  that  the  components  of  magnetic  field  that  are 

proportional  to  the  displacement  of  the  beam  do  not  diffuse  out  of  the  chamber  when 
2  2  2 

ojgtp  >>  1 (q  .  As  a  consequence  nsrfc  /a  does  not  increase  during  diffusion  and 

does  not  change  polarity.  Therefore,  the  drug  instability14  can  be  avoided  by 

choosing  the  various  parameters  to  give  w  >  o  at  the  commencement  of  the  diffusion 

3 


process. 


! 

During  the  main  acceleration  phase  the  significance  of  toroidal  effects  is 

reduced  because  v/y  ♦  o.  When  Ba  remains  approximately  constant  as  shown  in  Fig. 

o  “  f 

10,  or  increases  with  time  the  accelerated  beam  moves  closer  to  the  center  of  the 
minor  cross-section  of  the  torus.  This  may  be  seen  from  the  beam  equations,  which 
for  B  a  •  const. ,  are: 

O  o 

,  0  <6P  a> 

Ar  +  (y/y)  At  +  ui^Ar  -  (0oQ/y)  Az  +  — — ,  (40) 

Y  'mr0 

Az  +  (y/y)Az  +  id  ^ Az  -  -  (0  / Y)Ar.  (41) 

z  06 

With  the  exception  of  the  second  term  (y/Y)Ar,  Eqs.  (40)  and  (41)  are  identical  to 
Eqs.  (12)  and  (13).  For  n  *  1/2,  Eqs.  (40)  and  (41)  can  be  readily  solved. 

Introducing  a  new  variable  ■  Ar  +  iAz,  these  two  equations  can  be  combined  into  a 

single  equation 


iji  + 


l 

Y 


(;  +  ino9)  *  +  «o2 


♦ 


<6Pg> 


where  u 
o 


The  general  solution  of  Eq.  (42)  is 


(42) 


where 


,1/2  [(O^2  +4Uo2f« 
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Magn.  Field  (kG) 


and 


*i  / 


t  y  >mr  2 


w  (  ^ ^ ) 


+  I 


ft„fl  <5V 

(-^ - ^-)  *,  dt- 

Y  Ymr  1 
_ _ o _ 

w  (  ^  ,  ^2  > 


In  the  above  equations 


(0)  ^oQ  .  J,Q o9.2  .  2  .  ,  .  ,  , 

UJ+  -  J—  +  \(yr-)  +  u0  and  (Ktj,  <J>2) 


is  the  Wronskian 


2Y  ~  ,v2y 

of  the  two  Independent  solutions  ^  and  ^  of  Che  homogeneous  equation. 

Since  the  denominator  of  and  ^  increases,  when  the  betatron  field 
increases  in  time,  the  center  of  the  beam  moves  toward  the  center  of  the  minor 
cross-section  of  the  torus  during  the  acceleration. 


V.  Grad  B  Drift  in  the  Modified  Betatron 


Until  now  we  have  neglected  in  our  analysis  the  radial  gradient  in  the  toroidal 
magnetic  field,  BQ.  In  this  section  we  consider  carefully  the  effect  of  this  non- 
linearity  on  particle  motion,  assuming  that  it  is  the  dominant  non-linear  effect. 

The  linearized  equations  of  motion  of  a  particle  located  at  r  *  r^  + 

Ar  +  Sr,  z  *  Az  +  Sz  where  Ar  and  Az  are  the  beam  position  with  respect  to  the 
center  of  the  torus  (See  Figure  11)  are 


~  ,  r,  ^  r,  <$Pfi 

r.  +  (1  -  n)  ft  r.  -  n  ft  ( Sr  +  — tt-  Ar)  -  fta  ( 1 - -)  z,  +  ft  - , 

1  zo  1  s  zo  2  9o  r  1  zo  y  mr  ’ 

a  o  o  o 


r  r 

z.  +  nft  ^  z.  -  n  ft  ^(<5z  +  — ~  Az)  =  -  ft  (1 - -)  r  , 

lzolszo  2  9o  r  1 

a  o 


(43) 


where 


■  Ar  +  Sr, 
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Az  +  5z , 


z 


1 


n  -  n  /y  , 

zo  oz  o 


00 


Vv 


Choosing  n  *  j 
defining  ♦  ■  r  +  iz 


simplifies 

gives 


the  subsequent  analysis. 


Making  this  choice  and 


n 
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zo 
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r, 

(  S'!*  +  — — 
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'  '  lS„e  (1 


Re*! 
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+  n 


zo  Y  mr 
o  o 


(44) 


We  proceed  to  solve  (44)  perturbatively,  assuming  Reik  «  r  . 

1  o 

The  zero  order  equation,  neglecting  the  non  linearity,  is  easily  solved.  First 
an  average  is  performed  over  initial  positions  and  velocities  of  the  particles  to 
obtain  a  single  equation  governing  the  motion  of  the  beam  center,  At|>.  Denoting  this 
average  by  brackets  we  will  have  <5i|>>  ■  0  by  definition  and,  it  may  be  shown  that  as 
long  as  no  kinks  develop  in  the  beam  <6\p>  ■  0  ■  Once  the  averaged  equation  is 

obtained  it  may  be  subtracted  from  (44)  to  obtain  an  equation  for  5 .  Carrying  out 
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this  program,  we  find  for  the  zero  order  solution: 


Alp 


.(0) 


<5V 


~  ,1  rb  v 

nzo  ToMo  <2  "  T"  ns) 
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-i0>.  t  iltf-t 
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(45) 
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n  „  y  mr  (-r  -  n  ) 
zo  o  o  2  s 


-ioj.t  -iu>,t 

+  Ce  +  De 


where  A,  B,  C,  and  D  are  arbitrary  complex  constants  and  where  the  frequencies 
V..  are  given  by 

2  1 

2  ,  ,  ~  2,1  T_b  _  ^2 

2 
a 


"i,2-Tt5eoi(aeoi  +  45zoZ 

“  <3 


“3,4  ’I  i  (3to2  +  ‘V 


2  ,1 


We  shall  take  subscripts  1  and  3  to  correspond  to  the  +  signs. 

2 

I  r  b 

In  writing  (45)  we  have  assumed  that  neither  (y - y-  ng)  nor 

a 

(y  -  ng)  is  zero.  If  either  of  these  quantities  does  vanish  (corresponding 
physically  to  the  vanishing  of  net  radial  restoring  forces)  the  corresponding 
solutions  to  (44)  grow  secularly,  indicating  a  curvature  or  centrifugal  drift  in  the 
vertical  direction.  Below  we  assume  that  the  radial  restoring  forces  do  not  vanish 
for  either  beam  or  particle  motion.  In  addition  we  shall  make  the  assumption 
that  <5Pg>  *  0.  This  is  the  same  as  the  requirement  that  the  equilibrium  position 
of  the  beam  be  at  the  center  of  the  minor  cross  section  of  the  torus.  This 
assumption  has  no  effect  whatsoever  on  the  basic  physical  results  and  conclusions 
but  does  simplify  the  mathematics  somewhat. 


39 


We  may  return  to  (44)  and  calculate  the  first  order  correction  to  the  beam 
position,  The  equation  to  be  solved  is 


+  in.  di(1)  +  n  2 -  -5-  n  )  Ai|/  1 ) 


zo  2  2  s 

a 


Substituting  from  (45)  we  find  the  right  hand  side  of  (47)  is 


fn\  ( 2  2  -i(w.t  -  a.) 

<[Retp^  >  *  -i  I  I  “j  Pj  i^e  J  J  cos(a^t  -  c^) 


-?  J  “,«>]  >■ 
3-3  2  3 


where  we  have  defined 


A,  B,  C,  D  =  p  e  j  -  1,  2,  3,  4, 


and  where  we  have  assumed  that 


<C>  «  <D>  -  0. 


From  (47)  and  (48)  we  can  see  that  apart  from  oscillating  terms  the  net  effect 
of  the  radial  gradient  in  Bg  is  to  cause  an  outward  shift  in  the  equilibrium 
position  of  the  beam: 


V  4  9 


n  2  (i  _  b 

zo  ?  ~2“  ns ' 

a 


+  oscillating  terms. 
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This  result  is  simply  understood  as  the  result  of  a  balance  between  the  outward, 


"diamagnetic"  force  (which  tends  to  expel  the  beam  from  the  high  field  region)  and 
the  inward  radial  restoring  force. 

Since  the  values  in  (49)  depend  on  the  details  of  the  injection  process,  it 
is  difficult  to  draw  practical,  quantitative  conclusions  from  this  result.  However 
it  is  probably  safe  to  conclude  quite  generally  that  any  device  should  be  designed 

2 

1  r  b 

so  that  x  >>  — r—  n  for  all  times,  or 

2  2s 
a 


v/  Y  < <  4-  (  Y  Q  a/c)2, 

o  4  o  zo 


(50) 


for  the  case  of  a  perfectly  conducting  wall.  The  Yq  should  be  omitted  in  the 
parentheses  on  the  right  hand  side  of  (50)  in  the  case  of  a  poorly  conducting  wall 
(diffusion  time  short  compared  to  a  beam  oscillation  period).  Me  note  that  the 
constraint  in  (50)  is  independent  of  the  strength  of  the  toroidal  magnetic  field. 

Two  dimensional  computer  simulations  bear  out  our  claim  that  self  consistent 
beam  equilibria  exist  in  the  presence  of  a  gradient  in  Bg,  as  long  as  the  net  radial 
focusing  forces  (proportional,  basically,  to  the  denominator  in  (49)  but  generalized 
to  include  the  case  n  *  -j,  and  to  include  toroidal  corrections  to  the  self  fields) 
do  not  vanish.  In  Figure  12  we  show  a  succession  of  "snapshots"  of  a  beam  cross 
section,  which  remains  in  its  equilibrium  position  for  significant  times  compared  to 
r^/Vp,  where  VQ  is  the  single  particle  drift  velocity. 


D  9o 


p2/2v 


and  where  p  is  a  particle  gyradius.  No  drift  is  observed.  (This  is  not  just  a 
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Figure  12.  Four  snap-shots  of  the  electron  ring  minor  cross-section.  The 

values  of  the  various  parameters  used  in  this  computer  simulation 
run  are  shown  at  the  top  of  the  figure. 


,2 


visual  observation  but  is  obtained  from  a  plot  of  average  particle  position  vs. 
t ime . ) 

In  Figure  13,  on  the  other  hand,  we  present  a  case  in  which  net  radial  focusing 
does  nearly  vanish.  Loss  of  confinement  is  extremely  rapid  under  such  conditions. 
The  ring  drifts  vertically  with  an  average  speed  of  0.25  cm/nsec. 

VI.  The  Effect  of  Emittance 


Up  to  this  point,  we  have  dealt  with  the  equations  describing  the  motion  of  the 
center  of  an  electron  ring  in  a  modified  betatron  configuration.  In  this  section, 
we  discuss  the  effect  of  the  finite  emittance  on  the  equilibrium  of  the  gyrating 
electron  ring  using  the  beam  envelope  equation.  When  the  major  beam  radius  rQ  is 


I  V 

large,  n  *  y  ,  yq  *  constant,  —  «  1,  the  energy  mismatch  6yq  *  0  and  the  effect  of 

surrounding  walls  is  neglected,  the  beam  envelope  equation  in  the  paraxial 

13 

approximation  for  B  „  »  B  becomes  in  the  Larmor  frame  of  reference, 
o  o  oz 


+  “T  &T*  r-(s) 


oz 


2v/y„3  2 

B  2r  (s)  r  3(s) 

O  D  O 


(51) 


where  c  is  the  beam  emittance  (unnormalized),  s  *  0rQ  is  the  length  along  the  minor 

A 

axis  of  the  torus  and  r^  (s)  =  dr^/ds.  For  a  zero  emittance  beam  the  motion  of  the 
particles  is  laminar  and  the  equilibrium  is  called  either  laminar  flow  or  Brillouin 
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Figure  13.  Reducing  the  betatron  field  by  10  G,  but  keeping  the  remaining 
parameters  fixed,  as  in  the  run  of  Figure  12,  the  electron  ring 
drifts  rapidly  in  the  vertical  direction. 


flow.  The  equilibria  radius  in  this  case  is  obtained  from  Eq.  (51)  by  setting 


r 


b 


£  *  0 


and  is 


J. 

I  ~  2c_  ,2\\2 

“  no0  V  * 


(52) 


when  y  »  1.  For  a  finite  emittance  beam,  the  equilibrium  radius  can  be  determined 
o 

A  A 

from  Eq.  (51)  by  taking  r,  (s)  ■  0  and  is 
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Q  Y  2  —  — 
2c  .  C/ v  s2  ,  2(  o9  o>  i2i2 
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(53) 


For  small  envelope  oscillations  r,  »  r,  +  d,  with  d/r,  «  1  and  Eq.  (51)  gives 

b  b,eq  b,eq 


2[(no6/2V)2+  e2/tb‘eq]«.  0, 


that  has  a  period 


2r/c 


j*[«W*V>2  +  3  e2/0}  k 


(54) 


The  effect  of  emittance  on  the  equilibrium  of  the  ring  has  been  studied 
extensively  using  a  computer  simulation  code.  Numerical  results  from  the  computer 
simulation  are  given  in  Figs.  14  to  17.  For  a  beam  with  yq  »  7,  I  ■  10  XA, 
r  *  100  cm,  B  *  160G,  B  Q  *  1415G,  Eq.  (52)  predicts  that  the  equilibrium  radius 
for  the  Brilluin  flow  is  1  cm.  The  numerical  results  of  Fig.  14  give  also  a  radius 
of  1  cm,  that  for  all  practical  purposes  remains  constant  in  time.  In  this  run  the 
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Figure  14.  Root-mean  distance  of  the  electrons  from  the  center  of  the  ring  as 
a  function  of  time,  when  the  emittance,  e,  is  zero. 
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Figure  15.  Root-mean  distance  of  the  electron  from  the  center  of  the  ring  as 
a  function  of  time.  In  this  run  e  ■  50  mrad-cm. 
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Figure  16.  Four  3nap-shots  of  the  electron  ring  minor  cross-section  when 
50  mrad-cm. 
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Figure  17.  Root-mean  distance  of  the  elctrons  from  the  center  of  the  ring. 
The  envelope  oscillations  are  drastically  reduced  when  the 
toroidal  magnetic  field  is  increased  from  1415  G  to  1830  G. 


electron  beam  Is  injected  into  the  torus  with  a  rotational  frequency  that  is  half  of 

the  local  cyclotron  frequency.  Figure  15  shows  the  envelope  in  time  of  a  nor- 

rotating  beam.  The  various  parameters  in  this  run  have  the  same  values  as  those  in 

Fig.  14,  except  now  the  emittance  is  non-zero  in  the  Larmor  frame.  For  the 

equivalent  emittance  of  e  ■  50  mrad-cm,  Eq.  (53)  gives  an  r,  =*  1,21  cm  and  the 

b,eq 

numerical  results  1.24  cm.  Similarly,  Eq .  (54)  gives  period  of  2.13  nsec  and  the 

nunerical  results  2.28  usee.  Snap  shots  of  the  beam  cross-section  from  the  same  run 

are  given  in  Fig.  14.  The  oscillations  observed  in  the  run  of  Fig.  15  can  be 

avoided  by  "matching"  the  beam,  i.e.,  by  raising  the  B  Q  magnetic  field  to  1830G. 

*  oo 

This  value  of  magnetic  field  gives  an  r,  *1  cm,  which  is  the  radius  of  the  beam 

b,eq 

for  Brillouin  flow.  Numerical  results  from  this  run  are  shown  in  Fig.  17. 

The  electron  beams  discussed  so  far  in  this  section  were  monoenergetic  with 
finite  emittance.  Such  beams  have  an  axial  velocity  spread  equivalent  to  that  of  a 
cold  beam  with  energy  spread  Ay  that  is  given  by^ 


$  - 

b 


Actual  or  equivalent  energy  speed  in  the  direction  of  beam  propagation  has  an 
important  effect  on  the  dynamics  of  electrons  as  may  be  seen  as  follows.  The 
equations  describing  the  motion  of  individual  electrons  in  cylindrical  geometry  are 
identical  with  those  describing  the  motion  of  the  center  of  the  beam,  provided  that 

•v  2  -  2 

u>r  and  0)^  are  replaced  by 
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The  above  sec  of  equations  are  based  on  Che  assumption  that  toroidal  effects  can  be 
neglected.  In  addition,  it  should  be  emphasized  that  Ay  is  thermal  energy  spread 
and  not  the  energy  mismatch  6 yq  discussed  in  Section  III. 

For  n  *  y,  the  solution  of  the  individual  particle  equations  are 


(1  -  COSIDgt) 

r  -  r  -  6r  cosu)„t  -  5z  sin  u_ t  +  r - : - 

o  oB  o  B  oy  1 
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^  sinoigt 
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o  B  o  ^  o  y  1 
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(56) 


where 


n  b 

—  —■  (1/2  -  n  ). 

\  “e 


According  to  Eq.  (55),  when  ng  <<  1,  i.e.,  for  low  current  beams,  thermal  effects 

increase  substantially  the  minor  radius  of  the  beam.  In  such  beams  the  minor  radius 

varies  as  2r  Ay/y  • 
o  o 

In  contrast,  when  n  »  1,  i.e.,  for  high  current  beams,  thermal  effects  do  not 

S  ro  A 

change  significantly  the  minor  radius  of  the  beam,  which  varies  as  (— )  — . 

ns  Yo 

The  effect  of  the  axial  energy  spread  on  the  minor  cross-section  of  the  beam  in 
a  modified  betatron  geometry  has  been  studied  numerically.  Results  from  both  the 
high  and  low  current  beams  are  given  in  Fig.  18.  In  these  two  runs  the  various 
parameters  have  the  values  listed  in  Table  III.  In  the  high  current  beam  case  the 
minor  radius  expands  by  approximately  a  factor  of  two.  However,  in  the  lew  current 


51 


LOW  CURRENT  (100  A) 

A £/£  =  10% 

TIME  *  0.00  nsec 

rs 
r*i 


♦ 


N 

I 

I 


<N 

1  96.8  R  -►  10^2 


HIGH  CURRENT  (10  kA) 

A £/£  =  10% 


TIME  =  0.00  nsec 


1  96.8  103.2 


TIME  =  20.00  nsec 


TIME  =  10.00  nsec 


TIME  =  20.00  nsec 


Figure  18.  Snap-3hots  of  the  electron  ring  minor  cross-section  for  low  and 
high  current.  The  values  of  the  various  parameters  for  this  run 
are  listed  in  Table  II.  The  energy  spread  in  both  cases  is  10%. 


Table  III. 


! 


Parameters 


Low  Current 


High  Current 


0.1 


10 


Beam  Current  (kA) 

Beam  Energy  (MeV) 

Energy  Spread  (%) 

Major  Radius  (cm) 

Beam  Minor  Radius  (cm) 
Initial  Betatron  Field  (G) 
Toroidal  Field  (kG) 
External  Field  Index 
Self  Field  Index  (ns) 

Torus  Minor  Radius  (cm) 


10 

10 

100 

100 

1 

1 

116 

146 

1.4 

1.4 

0.5 

0.3S 

0.37 

23.4 

3.2 

3.2 

p 


M 
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beam  run  the  beam  expands  slgnif icantly  and  strikes  the  wall.  Therefore,  a 

substantial  energy  spread  can  be  tolerated  in  high  current  beams  without  a 

catastrophic  expansion  of  the  minor  radius  of  the  beam.  Such  energy  spread  may  be 

16—18 

required  to  stabilize  the  various  disruptive  instabilities. 


VII.  Summary 

In  this  paper  we  review  the  dynamics  of  ultra-high  current  electron  rings  in 

the  modified  betatron  configuration.  Our  discussion  addresses  mainly  the  evolution 

of  the  electron  ring  after  injection.  The  formation  of  the  ring  during  injection 

19 

has  been  analyzed  and  reported  previously. 

Our  work  includes  both  analytical  and  numerical  results  for  "cold"  and  "hot" 
rings.  The  conclusion  of  these  studies  is  that  equilibrium  states  of  ultra-high 
current  rings  in  a  modified  betatron  exist  over  a  wide  range  of  parameters.  These 
equilibria  are  realistic  and  accessible  with  state  of  the  art  injectors. 

The  results  presented  in  this  paper  are  based  on  several  simplified 
assumptions.  Among  them,  we  have  assumed  that  the  various  fields  are  free  of 
errors,  the  conducting  wall  that  surrounds  the  electron  ring  was  assumed  to  be 
smooth,  i.e.,  without  ports  and  gaps  and  the  external  field  index  was  assumed  to  be 
constant  in  time  and  space.  The  consequences  of  these  assumptions  is  presently 
under  investigation. 
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Appendix  A 

Diffusion  of  the  self  magnetic  field  of  the  beam  through  a  conducting 
liner. 

The  purpose  of  this  appendix  is  to  briefly  outline  the  calculation  of  the 
diffusion  of  the  self  magnetic  field  of  the  beam  through  a  conducting  liner. 
To  simplify  the  analysis,  it  is  assumed  that  the  electron  beam  is  located 
along  the  axis  of  a  straight,  circular  cylinder  of  inner  radius  a,  outer 
radius  b  and  thickness  6  »  b-a,  as  shown  in  Fig.  9.  In  addition,  it  is 
assumed  that  the  current  of  the  electron  beam  is  a  step  function  that  is 
turned-on  at  t=0.  Since  the  problem  of  interest  is  that  of  an  electron  ring 
inside  a  torus,  it  is  further  assumed  that  the  axial  inductive  electric  field 
goes  to  zero  at  r=rQ,  where  Tq  is  the  major  radius  of  the  torus. 

Neglecting  the  displacement  current  (quasi-static  approximation),  the 
fields  inside  the  cylindrical  conducting  shell  (a  <_  r  <_  b)  are  given  in  MKS 
units  by 

V  x  H  =  J ,  (A-l) 


V  •  B  -  0,  ( A-2) 

and 

J  *  oE,  (A- 4) 

+  +  + 
where  J  is  the  current  density,  H  is  the  magnetic  field,  B  is  the  magnetic 

induction  and  a  the  conductivity  of  the  conductor. 

.Assuming  that  all  quantities  are  independent  of  z  and  $ ,  Eqs.  (A-l)  to 

(A-4)  can  be  combined  to  give  a  diffusion  equation  for  the  azimuthal  component 

of  the  magnetic  field 
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if  [  f  If  (r  Vr.c>>]  -  .  (A-5) 


For  r  <  a  and  r  >  b,  the  conductivity  is  equal  to  zero  and  Eq .  (A-5)  becomes 


If  I  7  If  (r  H*  *  0  . 


(  A-6) 


with  the  only  acceptable  solutions 


(r>t)  ’  Iff  0(t)’  rb  <  r  <  a 


(A- 7) 


H*  (r,t)  -  H  x  (b,t)  b/r,  r  >  b 


(A- 8) 


where  6(t)  is  the  step  function,  and  is  the  beam  radius. 

To  complete  the  specification  of  the  problem  we  have  to  introduce 
boundary  conditions.  From  Eqs.  (A-3)  and  (A-8),  we  get 


3E  3H  (b ,  t)  , 

_ z  _  ..  d>  h 


which  a£ter  integration  yields 


3H^(b,t) 

Ez  (r,t)  -  y - ^ -  bin  —  , 

o 


(A- 9) 


where  we  have  assumed  that  E  (r,t)  is  zero  at  r»rn. 

Z  w 

Combining  Eqs.  (A-l),  (A-4)  and  (A-9)  at  r*b,we  obtain  the  first  boundary 
condition,  namely 


f T  (rH.)  1 .  -  y  (b , t)  b  in  . 


(A- 10) 


The  second  boundary  condi cion  is  furnished  from  Che  conCinuity  of  che  magnetic 
field  tc  r  ■  t,  i.e. , 


0(C)  -  HxCa.C). 


(  A— 11) 


Since  the  raagnecic  field  is  zero  aC  t-0,  Eqs.  (A-5) ,  (A-9)  and  (A-10)  take  the 
following  form  in  their  transformed  state 


i  ,  i 

72  +  7  sT  '  (  “2  +  ‘’"'’I  %  ’ 

dr  r 


[t  <-'.>]■ 


ovpH.(b.t)  b2in^~, 

41  ro 


( A- 1 2) 


(  A- 1 3) 


>(h ^  Cr.t)  I 
L  j  r*a 


(A-14) 


where,  the  Laplace  transform  of  is  defined  by 


H*  (r,p)  -  /  e~ptHd)  (r,t)  dt. 


Re  (p)  >  x. 


Equation  (A-12)  is  the  modified  Bessel  equation  and  its  solution  is 


H  ^  (r,p)  -  AIjCAr)  +  BKjUr)  , 


where  A  and  B  are  constants  and  X  *  opp. 

The  two  constant  coefficients  A  and  B  are  determined  from  Eqs.  (A-13)  and 


(A-14)  and  are: 


A-  (  2i^HVAb)  +  lb  • 

B-  f  I5I?4)Il0(Xb)  -  Xb  • 


where 


A-  Kq( Xb)  Ij  (Xa)+Io(Xb)Kiaa)+XbIn|-  [k :  (Xb)  Ij  (Xa) 


-  Kj(Xb)  IjCXb)]. 


(A-15) 


The  magnetic  field  in  the  region  a  £  r  <_  b  as  a  function  of  space  and  time  can 
be  obtained  by  inverting  the  Laplace  transform,  i.e. , 

c+i-  ePt 

Vr’° .  "  2^i  4-i-dP2^iA  {[K0(Xb)  +  Xb  Zn  ^Kl(Xb)J  ri(Xr) 

+  [lQ(  Xb)  -  Xb  in  Ix  ( Xb)  ]  Kx  ( Xr)  }  ,  (A-16) 


where  the  path  of  integration  is  a  vertical  line  in  the  complex  p-plane  to  the 
right  of  all  singularities  of  the  integrand. 

Equation  (A-16)  has  a  simple  pole  at  p»0  of  residue  a/r  and  an  infinite 
set  of  simple  poles  at  A-0. 

Contour  integration  of  Eq.  (A-16)  gives 


*  ik  jr  +  i . 

s-1 


\  2 

•  -*%  t/w0 
e 


Vasb,Jl(V>  *  W>  Vasb> 


*  «,b  [Ti(a.b,Ji<a.r)  -  Vv>  Ji<asb3]fA(if) 


cr*  a 


(A- 17) 
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where  u  ■  jJn,  A  ■  la,  Y,J  are  the  Bessel  functions,  a  are  the  roots  of 
u  s 

equation 


Y0(asb)Jl(asa) 


Yl(V)J0(asb)  +  °sb  *n  ^  tYl(asb)Jl(°sa) 


-  Yi(asa)J1(asb)]  -  0, 


(A-18) 


and 


s\da, 


a*a 


V  tY0(asb)J0(V°  ■  Y0(V)J0(asb)}‘  asb[Yl(asb)Jl(asa) 
-  VV)Jl(asb)]  +  °sb  £n  t  {baSfY0(txsb)Jl(°sa)  "  Yl(asa)J0(asb)^ 


+  V  tJ0(asa)Yl(asb)  '  Y0(asa)Jl(asb)^' 


For  a  a»  1,  Eq.  (A-17)  is  reduced  to 
s 


Vr-C)  ’  2l?  1  ,  e 

T  s=l 


-  a  t/ay 


^  [cos  a  (b-r) 


+  a  b  In  —  sin  a  (b 
s  rQ  s 


-*>]/(■ 


a  6  sin  (a  6)  -  a  b  in 

S  5  S 


—  cos  a  fij}, 

O  / 


(A-19) 


where  a  are  now  the  roots  of 
s 


cos  (a  6)  +  a  b  In  —  sin  (a  6)  *  0. 
s  s  r  s 

o 


(  A- 20) 


Equations  (A-19)  and  (A-20)  can  be  further  simplified  when  as5  <<  1,  In  this 
case 


59 


^(r.t)  "  2»«  ^  r  “  [co*  ag^b"r)  +  <*#b  to  sin  ag(b-r)] 


-a.  t/ayo  , 

*  }. 


which  at  r-b,  becomes 


-a  t/ op 

H^,(b,t)  -  yjjg-  {l  -  e  }, 


where 


2  1 


s  5b  tar. 


( A-21) 


Therfore,  to  lowest  order,  the  magnetic  field  in  the  region  r  >  b  is 


-  a  t/ oil. 

Vr’°  "  ‘HF  i1  “  •  h  r>b.  MKS  (A-22) 


or 


9T  ,  -a_2tc2/4TTO, 

H^Cr.t)  *  U  -  *  (•  CGS  ( A-23) 


Under  the  same  approximations,  the  electric  field  at  the  inner  edge  (r*a) 
of  the  conducting  shell  is 


Ez(a,t) 


-^°.2  «■&*) 


rft  -a  t/ou 


s  ww"0 


MKS  (A- 24) 


c  r 


y*.o  - 


I  2  .  ^O!.^.2''2'4"0 

"US  %  ta  (— J* 


CGS 
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Since  the  electric  field  is  unifora  in  the  region  r  <_  a,  Eq.  (A-24)  gives  also 
the  electric  field  that  acts  on  the  beam.  Substituting  Eq.  (A-24)  into  the 
energy  rate  equation 


me 


2  dy 
dt 


-e  v 


E, 


we  obtain  for  highly  relativistic  beams 


Ay 

y 


A-25) 


For  ^  *  10  and  *  0.1,  Eq.  (A-25)  gives  -^*0.46,  i.e.,  a  substantial 

reduction  in  the  energy  of  the  beam.  However,  for  highly  relativistic  beams 

v  ~  c  and  thus  the  current  of  the  beam  remains  approximately  constant, 
z 
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Appendix  B 


Description  of  the  Particle  in  Cell  Computer  Code 

The  NHL  modified  betatron  accelerator  is  designed  for  a  maximum  v/y  *  0.1, 
i.e. ,  the  current  is  high  enough  so  that  the  self-fields  of  the  ring  can  exceed  the 
externally  applied  fields.  A  realistic  theoretical  description  must  therefore  self- 
consistently  include  the  beam's  self  fields  as  well  as  the  effect  of  surrounding 
walls.  Since  this  is  difficult  analytically,  particularly  if  the  ring  is  displaced 
substantially  from  the  center  of  the  minor  cross-section  of  the  torus,  numerical 
simulations  are  useful  both  in  gaining  insight  into  the  important  physical  processes 
as  well  as  to  provide  a  method  to  check  the  applicability  of  specific  assumptions  in 
an  analytic  model. 

The  dynamics  of  the  accelerated  electron  ring  are  determined  by  forces  that 
vary  on  a  number  of  different  time  scales  which  range  from  the  electron  cyclotron 
period,  i.e.  ,  a  few  nanoseconds  to  the  beam  acceleration  time,  which  is  of  order  a 
millisecond.  The  code  described  here  is  taylored  to  simulate  efficiently  the 
various  phenomena  on  the  intermediate  time  scale.  This  time  scale  is  characteristic 
of  the  drift  (bounce)  motion  of  the  ring  after  equilibrium  has  been  established, 
rather  than  the  rapid  evolution  occuring  at  injection.  Simulation  of  a  single  turn 
around  the  major  axis  that  lasts  about  20  nsec  using  4K  particles  on  a  64  x  64  grid 
and  typically  takes  about  one  minute  on  the  NBL  Texas  Instruments  ASC. 

The  simulation  code  is  r-z,  spatially  two  dimensional,  i.e.  ■  0,  but  with 
three  velocity  components.  Although  Bg  is  used  in  calculating  the  particle 
trajectories,  it  is  not  solved  self-consistently ,  i.e.,  is  assumed  to  be  generated 
from  external  coils  only.  This  assumption  is  valid  to  first  order  in  v/y.  The 
radiative  term  (displacement  current)  is  also  ignored,  i.e.,  the  code  uses  the 
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Darwin  model  for  Maxwell' s  equations 


The  electrostatic  potential  is  computed  from  Poisson's  equation 


V2 


p/  e 


(B-l) 


and  the  magnetic  vector  potential  from 

’V7--W  (B-2) 

with  the  boundary  condition  $  ■  Ag  ■  0  at  the  conducting  wall. 

Equations  (B-l)  and  (B-2)  are  solved  by  Fourier  decomposition  in  the  z- 

direction  and  then  by  Gaussian  elimination  of  the  resultant  tridiagonal  matrix  of 

2 

equations  obtained  from  a  3-point  differencing  scheme  for  .  The  inverse  Fourier 
transform  yields  AQ  and  *  on  the  grid.  Note  that  the  B  particle  velocities  are 
advanced  using  the  conservation  of  cannonical  momentum  in  the  9  direction,  the 
equation  for  Ag  is  therefore  not  properly  time  centered  since  the  velocities  from 
the  previous  time  step  are  used  to  calculate  the  currents  from  the  cannonical 
momenta.  This  method  was  chosen  primarily  for  its  speed  and  simplicity  but  care 
must  be  taken  in  applying  the  code  when  the  inductive  acceleration  of  particles  in 
the  9  direction  is  significant. 

If  boundary  conditions  other  than  Ag  or  $  equal  zero  on  a  rectraqular  grid  are 
desired,  it  is  possible  to  obtain  relatively  arbitrary  boundary  conditions  using  the 
capacitive  matrlx(Buneman)  technique.  In  this  method  a  matrix  is  generated  which  is 
the  Green's  function  for  discrete  "wall"  points  within  the  system.  Then  at  each 
time  step  the  field  solver  described  above  is  used.  The  potential  at  the  discrete 
(wall)  points  is  obtained.  By  multiplying  this  vector  potential  by  the  invense 
matrix  obtained  previously  a  set  of  (wall)  currents  or  changes  is  generated,  which 
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is  used  co  specify  che  wall  potential.  These  "wall''  sources  are  then  added  to  the 
original  beaa  source  and  the  field  solver  is  used  again.  The  result  by  the 
principle  of  superposition  is  correct  inside  the  system  and  has  the  correct  boundary 
condition  on  the  "wall". 

Particle  Push 

The  motion  of  the  electrons  is  governed  by  the  Lorentz  force 


d(  yV) 
dt 


^  (E  +  V  x  B), 


>  «► 

where  E  and  B  are  the  total  electric  and  magnetic  field  respectively. 

In  component  form  the  equations  used  to  update  the  velocities  and  positions  at 
each  time  step  are 


6 


d(Wr) 

dt 


(E  +  V  x  B  -  V  x  Bj 
m  r  0  z  z  0 


d(W  ) 
2 

dt 


_! 

m 

o 


<Er 


dr 

dt 


dz 

dt 


-  V  V  V  Be) 


-  V 

r 


-  V 

z 


(B-3) 


CB-4) 


(B-5) 


(B-6) 


CB-7) 


Equation  (B-3)  is  used  at  each  timestep  to  compute  y  Vg.  Equations  (B-4)  and  (B-5) 
are  coupled.  To  advance  these  velocities  a  leap-frog  scheme  is  employed.  At  time 
t  all  fields,  positions  and  Vg  are  known  exactly.  Vf  and  V2  however  are  known  at 
t+1  and  the  velocities  will  be  advanced  to  t+l/2.  Before  the  equations  are 
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differenced,  it  is  convenient  to  rewrite  the  equations  in  terms  of  the  relativistic 
momenta.  Letting 


U  -  YV 


the  equations  (B-3),  (B-4)  and  (B-5)  can  be  written  as 


V  TT  ~i~  A0 

O  0 


dur  a  U  L 

IF  “  it  (^r+  V  V  V  B0)  +  IT 


dU 


(YEz  "  V  Br  +  V  V 


dr 

dt  "  Y 


.  u 

dz  z 

dt  *  y 


This  formulation  is  then  differenced  by  substituting 


U  +  1/2  - 
dt 


for  dt 


and 


♦t+l/2  +  tft-1/2 


for  U  . 


(B-3a) 

(B-4a) 

(B-5a) 

(B-6a) 

(B-7a) 


After  making  these  substitutions  in  equations  (B-4a)  and  (B-5a)  it  is  straight¬ 
forward,  if  somewhat  tedious,  to  solve  the  two  coupled  equations  for  Uzc+^^  and 
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Urc  1  .  Since  Ur  and  Ut  are  calculated  at  the  half-times teps  while  Ug  is  known  at 
the  full  timesteps,  Y  is  not  known  at  time  t.  This  difficulty  is  overcome  using 
iteration. 

After  advancing  the  velocities  the  particle  positions  are  advanced  using  a 
simple  centered  difference 


r 


t+1 


9 


and 


t+1  _  t  x_  Uz 
s  *  z  +  5T .  ^  « 


Once  all  the  velocities  and  positions  are  found  the  new  current  and  charge  densities 
are  used  to  update  the  fields  and  a  new  time  step  begins. 
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